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Abstract 

The controllability of charge transport through an organic molecular spin-valve 
system is theoretically investigated on the basis of a Su-Schrieffer-Heeger model 
combined with the non-equilibrium Green's function formalism. We show how the 
formation of polaron in the organic sub-structure leads to a hysteretic conductance 
switching, via sweeping either the bias voltage or the electrochemical potential. We 
further obtain an exponential dependence of the magnetoresistance as a function 
of the applied bias voltage. The implications of calculated results in relation to 
experiments and device applications are addressed and commented. 
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1 Introduction 



Due to their compatibility in processing and tunability in electronic proper- 
ties, 7r-conjugated organic semiconductors (OSEs) have been utilized broadly 
for various device applications [1,2,3,4,5]. Recent progress in the field of OSE 
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electronics has revealed a number of technologically important transport char- 
acteristics, such as current rectification [6], Coulomb blockade [7], Kondo res- 
onance [8], negative differential resistance (NDR) [9,10,11,12], bistable switch- 
ing [13,14,15,16,17,18,19], and giant magnetoresistance [20,21]. 

An OSE system is characterized by not only its electronic structure but also its 
strong electron-phonon (e-ph) coupling that induces polarons [22,23]. In view 
of charging-induced conformational change [14,24], the formation of nonequi- 
librium polaron should be a reasonable mechanism of the conductance switch- 
ing and/or NDR in OSE electronics. This idea is first shown theoretically 
by Galperin et al. [25] with a simplified model, in which the molecular wire 
consists of a single electronic state coupled with a single Brownian oscillator. 

Apparently, the Su-Schrieffer-Heeger (SSH) model [26] renders a more realis- 
tic description of an OSE system. It captures the essential characteristic of a 
conjugated molecule, where the strong e-ph coupling leads it to the polaron 
charged states and dimerized ground state. Indeed, the SSH model has been 
remarkably successful in the study of the electronic conductivity and opti- 
cal phenomena in 7r-conjugated OSEs. In addition to conducting polymers, 
the SSH model has also been extended to molecular systems [27,28], carbon 
nanotubes [29], and DNA molecules [30]. We have recently [31] exploited this 
model for an OSE system, and demonstrated that a NDR may appear due 
to the annihilation of a predoped polaron; i.e., an equilibrium charge transfer 
state that may exist before applying the external bias voltage. The resulting 
NDR can also be spin-polarization sensitive, leading to a magnetoresistance 
significantly larger than 100%. 

In this work, we elucidate the possibility of a polaron mechanism for the 
bistable switching and hysteresis in OSE nanostructures, and further for the 
giant magnetoresistance in OSE spin- valve devices where ferromagnetic metal 
electrodes are used. More importantly, we highlight that these technologically 
important properties are likely achievable by a wide range of OSE systems. To 
do that, moderate/typical values of the SSH parameters will be adopted in this 
work to avoid the formation of predoped polarons mentioned earlier. Besides 
the model and its parameters for the OSE spin-valve nanostructure of study, 
Sec. 2 presents also the transport theory based on the Keldysh nonequilibrium 
Green's function (NEGF) formalism [32,33,34], together with its implementa- 
tion to the evaluation of current in response to linear sweeping bias voltage. 
Hysteretic conductance switchings controlled via bias voltage and via an elec- 
trochemical means are then numerically investigated in Sec. 3 and Sec. 4, re- 
spectively. The giant magnetoresistance property of OSE spintronic structures 
is addressed in Sec. 5. Finally, we conclude in Sec. 6. 
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2 Model and theory 



2.1 A Su-Schrieffer-Heeger static polaron model 

We adopt a non-degenerate SSH Hamiltonian for the isolate OSE system 
[20,30], 

#OSE = J2{ e ° C n,v Cn >° ~ ( to ~ ( _1 ) nt l ~~ a o(«n+l ~ U n)] 
n,a 

X «a C n+l,- + C n+l,* C n,*)\ + ~^ £(«n+l ~ (1) 

Here, c+ CT (c„ jCr ) denotes the creation (annihilation) operator of an electron at 
the n th site with spin a of the OSE system. The e-ph coupling induced on-site 
lattice displacement parameters {u n } are treated classically and will be deter- 
mined in a self-consistent manner to be described later (cf. Sec. 2.2). The input 
parameters involved in Eq. (1) are the on-site energy e a , zero-displacement hop- 
ping integral t , non-degeneracy parameter ti, electron-phonon coupling a D , 
and lattice elastic constant K . 

For the electrodes, we choose two symmetrical 3d ferromagnetic (FM) tran- 
sition metals; each of them is described by a one-dimensional tight-binding 
Hamiltonian with the Stoner-like exchange term [35], 

H FM = ^2\ e fdn,adn,a + tf{d^ a d n+lt(7 + d+ +1 ^d nt0 .) } Jfi^n^n,] ~~ ^n^n.-l) ■ 

(2) 

Here, d+ (d n>a ) denotes the creation (annihilation) operator of a 3<i-electron 
at the n th site with spin a of the metal electrode; e/ is the on-site energy of a 
metal atom, tf the nearest neighbor transfer integral, and J/ the Stoner-like 
exchange integral. 

The coupling between the OSE and each of the FM electrodes is described by 
a tight-binding hopping integral, 

t c = /3(t f + t ). (3) 

Here, (3 is the OSE-FM binding parameter, which for simplicity is set to be 
spin-independent. However, charge transport can be spin dependent. Accord- 
ing to the two-current description, the total current arises from the majority- 
spin (up-spin) and the minority-spin (down-spin) of the 3d electrons of the 
FM electrodes [36]. We also neglect the spin-flip processes during the spin- 
dependent charge transport that takes place between the FM electrodes. The 
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above simple model for the molecule-metal coupling allows us to focus on the 
effects of OSE molecular system on the spin-dependent transport. 



2.2 Non- equilibrium Green's function formalism 



The quantum transport, in terms of the conductance switching, hysteresis and 
giant magnetoresistance behaviors of the model FM/OSE/FM system, will be 
calculated by means of the NEGF approach, based on the Keldysh formalism 
[5,32,33,34]. To proceed, the infinite FM/OSE/FM system is divided into three 
distinct regions: The central scattering region (S) that consists of the OSE 
together with a few metal atoms attached to each of its ends [34], and two 
semi-infinite FM electrodes (L and R) that serve as charge reservoirs and also 
set the temperature and electron distribution of the junction in the steady 
state. As results, the Hamiltonian for the whole system is mapped onto a 
3x3 block matrix, 



H 



total 



Ht. H 



LS 







Hsl H s Hi 



SR. 



V 



H KS H R 



(4) 



with i?L, Hs, and Hr being the isolated Hamiltonians for the L-electrode, the 
central scattering S-region, and the R-electrode, respectively, while i^s — H$ L 
and i^Rs = i?sR the couplings between them. 

The applied electric potential bias V does not change the electronic structure 
of L and R leads, but just shifts their energy levels by eV/2 and — eV/2, respec- 
tively. This leads to the electron distribution function f(E, fi a ) = [e ( - E ~^ a ^ kBT + 
where /xl/r = /x eq ± eV/2, with /x eq = Ep being the Fermi energy of the 
FM electrode. In contrast to the electrodes, the effect of applied bias on the 
S-region is characterized by the electric potential function <&(x). With the 
boundary conditions of $(xi) = V/2 and &(x Ns ) = —V/2, it shifts the n th - 
site energy by — e$(x n ) in the S-region. As it is correlated with the electronic 
and lattice structures of OSE, <§>(x) should be evaluated together with the S- 
region density matrix and lattice distortion in a self-consistent manner, which 
will be elaborated soon [cf. Eqs. (10)— (12)]. 

Eliminating the L- and R-electrodes degrees of freedom results in the NEGF 
of the central S-region [33], 



G(E) = 



E-H s + e$(x) - £ L (£) - £ R (£) + *0 H 



(5) 
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The effect of L- or R-electrode coupling is incorporated via the self-energy, 

E a (E) = H Sa - 1 H aS ; a = L, R. (6) 

In particular, T^/ R (E) = z[El/r(£') — T, R / L (E)] amounts to the electrode 
coupling-induced broadening function. 

The electron transmission probability function reads in terms of the NEGF at 
a finite bias as [33] 

T(E, V) = Tr[r L (E)G(E)r R (E)G\E)\. (7) 
The electric current can then be evaluated as [33,37,38] 



2e 
~h 



oo 

J T(E,V)[f(E,p L )-f(E,p R )}dE. (8) 



We shall show later that the above expression remains practically valid in the 
evaluation of current under the quasi-steady-state voltage sweeping condition; 
cf. the comments to the end of this subsection. 

The NEGF determines also the density of states D(E) and the density matrix 
p in the S-region [34]. They are given respectively by 

D{E) = ^\G{E) - G\E)i (9) 



and 



oo 

P=^ E / f{E^ a )[G{E)T a {E)G\E)]dE. (10) 



The number of carrier electrons in the S-region is iV = Tr (pS), with S denot- 
ing the overlap integral matrix that reduces to the unit matrix if an orthogonal 
basis set is used. 

We are now in the position to complete the self-consistent evaluation of the 
NEGF. The involving electric potential function $(x) in Eq. (5) is related to 
the charge density, and the latter depends further on the steady-state lattice 
distortion at a given finite applied electric potential bias. 

Let us start with the Poisson equation that relates $(x) and the on-site charge 

\Pnn) i 

V 2 ^(x) = -J2—S(x n -x). (11) 

n E 

Here, £~o denotes the vacuum permittivity, and x n = (n — l)a + u n , with a 
being the e-ph coupling-free lattice constant and set to be 0.122 nm. The e-ph 
coupling-induced on-site lattice displacement is of the boundary conditions of 



5 



u\ = un =0. To solve the Poisson equation we introduce = (0.5 — 

x/L)V + 5<&(x), with L denoting the S-region length. 5Q(x) satisfies the same 
Poisson equation, but with the boundary conditions of 5$(xi) = 8$(xn o ) = 0. 
It can be readily evaluated in terms of the on-site charge {p nn } via the method 
of images, in which 5Q(x) is a sum of the individual image contributions from 
all point charges. To avoid the infinity associated with point charges, we follow 
the Appendix 3 of Ref. [39] and analytically integrate the usual electrostatic 
potential within the small area around each point charge separately. 

Consider now the lattice distortion, which plays an important role in de- 
termining the properties of OSE due to the strong e-ph interactions. For a 
steady or quasi-steady state, the lattice distortion is assumed to catch up 
with the variation of charges under a slow voltage sweep condition treated in 
this work. In this case, the OSE lattice distortion can be determined in terms 
of the density matrix elements via the Hellman-Feynman variation theorem of 
= dEs/dun = dTr(H s p)/du n . It is that 

2a (Pn,n-l - Pn+l,n) + K (2u n - M„_i - U n+1 ) = 0. (12) 

Equations (10)-(12) constitute a set of coupled equations for solving p, 
and {u n } in a self-consistent manner by iteration. The converged result in the 
electric potential $(x) is then used to determine the NEGF of Eq. (5), and 
then the I-V characteristics of Eq. (8). The most time consuming part of the 
aforementioned self-consistent evaluation is Eq. (10). The involving integrand 
is non-zero over a broad energy range and often has sharp structures that 
makes Eq. (10) numerically challenging. To facilitate this problem, let 

P = Pcq + 5p, (13) 

where 

oo 

p eq = - E / f(E,p cq )[G cq (E)T a (E)Gl q (E)]dE, (14) 

«=L,R„ 00 

is the equilibrium density matrix, and Sp the non-equilibrium correction. The 
evaluation of p eq requires a broad energy range and couples only with Eq. (12) 
as V = is implied. The involving integration will be tackled by taking the 
advantage of the contour integration technique [34,40]. On the other hand, 
the non-equilibrium correction 5p involves the self-consistent determination, 
together with Eqs. (11) and (12). It will be evaluated directly, as the integra- 
tion is now required only for a limited energy range, where the integrands in 
Eqs. (10) and (14) differ significantly. This is the same situation as the evalu- 
ation of current [Eq. (8)], and can thus be carried out with a fine grid. Note 
that p cq is naturally real but 5p may not, unless a real basis set functions is 
employed. This is the case of this work, and the resulting p does assume a real 
matrix. 
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To conclude this subsection, let us comment on the quasi-steady-state condi- 
tion that will be applied in this work for the evaluation of I-V characteristics 
subject to a sweeping voltage of V(t) — V ± V't. Here, V = or V max for the 
sweeping up (+) or down (— ) direction, respectively; V > denotes the con- 
stant rate of bias sweeping, which is typically 0.01 ~ 0.1 V/sec in experiments. 
In our calculations of the quasi-steady-state current of interest, we set V = 0.1 
V/sec. Note that all experimental observations of the correlated electron and 
lattice degrees of freedom on the OSE are functionals of density matrix p. As 
the bias sweeping is linear, one can map the time evolution of density matrix 
onto the bias sweeping one, p(t) — > p(V). Starting from p(t = 0), the density 
matrix evolution is then evaluated via p(t + At) = p(t) + (Ap(t))r, where the 
second term denotes the time average of Ap(t) over the range of [t,t + T\. 
with T being the measurement duration (set to be T = 1 second) under the 
quasi-steady-state condition. The corresponding I(t) — > I(V) in each of the 
measurement time interval is then formally the same as Eq. (8). 



2. 3 Model parameters 

Reasonable values for the model parameters must be determined from other 
sources, such as the band structure calculations and/or experimental current- 
voltage (/-V) characteristics. In this study, the FM electrodes are chosen to 
be cobalt metal. With its experimental 3<i-band structure and Fermi energy 
of E F = —4.9 eV [21,41], the tight-binding parameters for the Co metal in 
Eq. (2) are determined as tf — 1.5 eV, J/ = 1.45 eV, and e/ = —6.5 eV. The 
parameters of the OSE in Eq. (1) are chosen as e Q = —5.5 eV, t = 2.5 eV, 
fi = 0.04 eV, a Q = 5.0 eV/A, and K = 21.0 eV/A 2 . These are the modified 
values from those of a typical conducting polymer, but lead to a relatively large 
band gap, which is about 2.8 eV for the OSE length of N Q = 10 sites considered 
in this work. In a realistic OSE molecule, the on-site energy e G can be site- 
dependent and also can be adjusted by applying a gate voltage [4]. The OSE- 
Co metal binding parameter (3 affects the magnitude of current sensitively, but 
it is hardly controlled due to the weak interaction. In the calculations, one can 
adjust it so that the calculated I-V curve is of about the same amplitude of the 
measured one. We find that choosing (3 = 0.25 for the conductance switching 
via applied voltage sweep and (5 = 0.5 for that via electrochemical potential 
sweep will produce the current in the magnitude of nA for some nano-size 
molecules reported experimentally [13,14,15,16,17,18,19]. 

For a self-consistent evaluation of quantum transport based on the NEGF for- 
malism, the central scattering S-region is chosen such that its Green's function 
can properly represent the bulk property at the boundary. For the coupling pa- 
rameters given above, it is sufficient to have the S-region a total of 20 + 10 + 20 
sites; i.e., to include 20 Co electrode atoms to each end of the model OSE of 
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Fig. 1. Density of states of Co 3d band and OSE structure at 11 K. The solid 
(dashed) lines correspond to majority (minority) electronic bands. The thick (thin) 
curves are for OSE (Co), with the phenomenological Lorentzian width of 0.15 eV 
being included. 

length N Q = 10. For the OSE-Co binding parameter used in this work, it is 
found that the OSE subsegment is charge neutral in the absence of applied 
bias (V = 0). Changing the relative energy levels of OSE with respect to 
electrodes may lead to intra-regional charge transfer, forming a preexisting 
(charged) polaron state in the OSE before applying bias voltage [31]. 

Figure 1 depicts the band structures of the OSE and the Co metal electrodes at 
a low temperature (T = 11 K). This figure summarizes the model parameters 
adopted in this work (except for (3 which in this figure is set to zero). For 
the OSE, the evaluated band gap is about 2.8 eV, which well reproduces the 
band diagram of the organic spin- valve system in a recent experiment [21]. 
For the Co electrodes, our calculation is consistent with the real material [41], 
showing a completely filled 3d band and a half-filled 3d band for the majority- 
spin (up-spin) and the minority-spin (down-spin) electrons, respectively. 

Quantitative simulations of experiment results [21,42,43,44] may need to in- 
corporate the spin dynamics, such as spin-flip, diffusion and accumulation, 
which are however not included in this work. Despite of its simplicity, we shall 
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show soon that the present model possesses some interesting properties that 
can highlight the important transport mechanism in relation to experiments. 



3 Hysteretic conductance switching via bias voltage sweep 

In this section we report the hysteretic bistable switching I-V behavior of the 
model Co/OSE/Co structure. In contact with experimental reality [13,18], the 
bias sweep rate in our calculations is fixed at a representing value of 0.1 V/sec, 
at which the quasi-steady-state condition is reasonably satisfied. 

There are two distinct transport measurement configurations, parallel (P) and 
antiparallel (AP), with respect to the relative magnetization orientation of 
two FM electrodes. Generally, in the two-current model, both the majority- 
majority and minority-minority (or majority-minority and minority-majority) 
transports are permitted in the P (or AP) configuration. However, for the 
cobalt electrodes in the present work, the full-filled up-spin electrons can- 
not transport in P-configuration under a positive bias voltage V > and 
the charge carriers are only the half-filled down-spin electrons. Contrary, the 
charge carriers in AP-configuration under V > are only the up-spin electrons, 
driven from the full-filled majority-spin band of L-electrode to the half-filled 
minority-spin band of R-electrode of opposite magnetization orientation. 

Figure 2 depicts the resulting hysteretic I-V characteristic in P-configuration. 
The bias sweep directions are indicated with the arrows. The OSE nanos- 
tructure is shown to have an insulating state and two conducting states in 
the indicated sweeping bias range. These conducting states are switched on 
(off), respectively, at 1.5 ± 0.2 V (0.75 ± 0.2 V) and 2.6 ± 0.1 V (2.1 ± 0.1 V) 
in the sweep up (down) direction. The voltage-triggered two-step hysteretic 
conductance switching has been observed experimentally in some molecular 
nanojunction systems [18], often associating with voltage-triggered molecular 
confirmation change. On the other hand, that observed in Fig. 2, according 
to the e-ph coupling Hamiltonian used here (cf. Sec. 2), can only be of po- 
laron in nature. The insert in Fig. 2 shows the sweeping bias-induced charge 
modification Aq D within the OSE sub-structure. The synchronization between 
the Aq Q -V and I-V curves demonstrates clearly that the two voltage-triggered 
conducting states have 1.15e and 1.44e charges, respectively. These are neg- 
atively charged polaron states, with the clear evidence of the e-ph coupling 
induced lattice distortion being presented soon (cf. Fig. 3). 

To further elaborate what happens to the conductance bistability of the OSE 
nanostructure, Fig. 3 depicts the transmission function and lattice distortion 
for the bistable states, exemplified with the P-configuration at V = 1 eV (cf. 
Fig. 2). The resulting window between /x L and /xr for the electric transmission 
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Voltage(V) 



Fig. 2. The hysteretic current as the function of sweeping bias potential for the 
model OSE in junction with Co electrodes in the parallel configuration. The linear 
voltage sweep rate is at 0.1 V/sec. The insert is the bias voltage-induced charge 
variation in the OSE substructure. The OSE-Co metal binding parameter (3 = 0.25. 

T(E), as indicated in Fig. 3a, clearly distinguishes between the conducting 
state (solid curve) and the insulation state (dash curve). Note that the lattice 
distortions in Fig. 3b are reported in terms of the e-ph coupling induced bond 
length variation y n = u n+ i — u n between nearest neighbors. The dashed curve 
there for the insulating state is effectively the same as the equilibrium (V = 0) 
lattice distortion. In contrast, the solid curve in Fig. 3 demonstrates that a 
charged polaron state is formed in the conducting state toward the right-side 
of the OSE-substructure due to the external bias-induced spatial symmetry 
broken. 

The polaron mechanism for the observed hysteretic conduction switching be- 
havior in Fig. 2 can now be summarized as follows. In the equilibrium state 
(V = 0), the OSE substructure is charge neutral, with its lowest-unoccupied 
and highest-occupied molecular orbital (LUMO and HOMO) electronic wave 
functions delocalized over the entire S-region. The chemical potentials of elec- 
trodes are the same as the Fermi energy, /il = //r = Ep, locating in between 
the LUMO and HOMO levels of the OSE; no current is observed and the 
OSE is in an insulating state. Applying a positive bias V > 0, which elevates 
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Fig. 3. Bistable switched-on (solid-curves) and switched-off (dashed-curves) states 
for the P-configuration at bias voltage V = 1.0 V. (a) Transmission coefficient of 
down spin electrons. The vertical dot lines indicate the window between [1l and 
fiR for the charge transmission; (b) Lattice distortion y n = u n+ \ — u n of the OSE 
substructure. 

/il = E F + eV/2 and depresses = E F — eV/2, drives electron transport from 
L-electrode to R-electrode through the unoccupied levels of OSE substructure. 
Since the two bistable conductance states in Fig. 2 are of the same physical 
origin, let us focus on the first one for illustration. As the sweep-up voltage 
approaches to the vicinity of the LUMO level, electron starts to migrate into 
the OSE, followed by the e-ph coupling-induced reorganization that leads to 
the formation of polaron. In other words, V on = 2(_E LUMO — E F )/e defines 
the switch-on voltage of the first conducting state. However, as the OSE gets 
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Fig. 4. Same as Fig. 2 but for the antiparallel configuration of the Co electrodes. 



charged, the e-ph coupling-induced reorganization process occurs, forming a 
polaronic gap state that serves as the first conducting state in Fig. 2. The fact 
that -Epoiaron < -^lumo due to the reorganization accounts therefore for the ob- 
served hysteresis, since V q q = 2(_E po i aron — Ep)/e is the switch-off voltage in the 
down-sweep direction. In the case of study here, V on as 1.5 V and V Q g pa 0.75 V 
(cf. Fig. 2). The reorganization energy of the negatively charged polaron gap 
state can be obtained as £"lumo — -^polaron = e(K n — V Q s)/2 pa 0.37 eV. 

Figure 4 depicts the hysteretic I-V characteristics and the corresponding 
charge modifications for the AP-configuration of the electrode magnetization 
directions. Similar to the P-configuration, the two-step hysteretic conductance 
switching behavior appears also here. This observation indicates that the hys- 
teresis itself is originated from the OSE due to the polaron formation, re- 
gardless the contact is magnetic and nonmagnetic in this aspect. However, 
in the present Co/OSE/Co model study, the hysteresis curves in the AP- 
configuration are narrower than those in the P-configuration. This difference 
gives us an additional handle by using a magnetic contact on the conductance 
control since one can alter from the P- to the AP-configuration of magnetiza- 
tion by applying an external magnetic field, especially when the electrode is 
made of soft magnetic material. 
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4 Hysteretic conductance switching via electrochemical potential 
sweep 

We now turn to another type of conductance switching, controlled by sweeping 
the Fermi energy [or e/ effectively in eq (2)] of both electrodes simultaneously 
via an electrochemical means, or equivalently sweeping the gate voltage via 
an microelectronic means that adjusts the on-site energy of the OSE. Here, 
the bias voltage is kept fixed at a small value to open a transmission window 
between /i L and /xr. For instance, in the experiment of He et al. [18], the cur- 
rent through Au-polyaniline-Au nanoelectronic junctions was measured under 
a linear sweep of the electrochemical potential at a fixed bias voltage of 20 
mV. The fast conductance switching behaviors of the electrochemical system 
is closely related to chemical and biosensor applications [18,45], while the gate 
voltage modulation is of great importance for the organic transistors and mem- 
ories [16,46]. They share the same mechanism; i.e. the control of the polaron 
(reduced/oxidized) state of the molecule by electrochemical potential or gate 
voltage, and can therefore be qualitatively illustrated with the same model. 
As discussed earlier, the polaron formation is responsible to the conductance 
switching. It arises from the coupling between the electronic states of the OSE 
near the LUMO energy and that of the metal leads near (j, a (a = L, R). When 
li a approaches to the polaron energy level (which is lower than the LUMO of 
the ground state) an abrupt change of current occurs. 

In the following calculations, we fix the electrochemical potential sweeping rate 
at 0.1 eV/sec that supports the quasi-steady-state condition. As described to 
the end of Sec. 2.2, Eq. (8) together with Eqs. (10)-(12) are used again for the 
evaluation of the quasi-steady-state current as the time-averaged one in the 
interval of every second. 

Figure 5 shows the dependence of the charge transport current on a sweeping 
of the electrochemical potential Ep, in terms of its difference from the initial 
Fermi level Ep (= 4.9 eV) of the metal leads in the P-configuration of elec- 
trodes. The bias voltage between the electrodes is kept at a constant value 
of 20 mV. The OSE-Co metal binding parameter is /3 = 0.25 in Fig. 5a and 
P = 0.5 in Fig. 5b, respectively. While the value of (3 affects the magnitude of 
current, the hysteretic behavior of conductance switching is shown to be an 
intrinsic character of the OSE, regardless the details of OSE-electrode cou- 
pling. Making in contact with the experimental output of current, which is of 
nA in magnitude for a nano-sized molecule [16,18], (3 = 0.5 is favored for the 
calculations of the potential-dependent conductance switching behavior. 

As shown in Fig. 5b, enhancing the potential leads to a rapid increase of the 
current during Ep — Ep ~ 0.5 eV and 0.75 eV, as the OSE switches from 
the insulating to the conducting state. The hysteresis behavior here is at- 
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Fig. 5. The hysteretic current as the function of sweeping electrochemical potential 
for the model OSE in junction with Co electrodes in the parallel configuration, at 
the fixed bias voltage of V = 20 mV. The OSE-Co binding parameter (a) (3 = 0.25 
and (b) f3 = 0.5, respectively. The initial Fermi level Ef = 4.9 eV. 

tributed to the formation of a bipolaron state that localizes two electrons in the 
OSE region. The formation of a bipolaron (not depicted here) is further con- 
firmed with various OSE chain lengths of iV = 10 ~ 40. The electrochemical 
potential-dependent charge transport behavior is similar to its bias voltage- 
dependent counterpart. In both cases, the same physical mechanism; i.e., the 
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polaron-assisted charge transport is operational. A similar phenomenon is also 
found for the AP-configuration of electrodes. 

The potential dependent spin-polarized charge transport through OSE has 
been studied theoretically by Xie et al. [35] in the context of the ground-state 
properties of ferromagnetic metal/conjugated polymer interface. The conclu- 
sions there are that electrons can transfer to the polymer and form bipolarons 
from the magnetic layer through an interfacial coupling by adjusting the rel- 
ative chemical potentials of the contact and the polymer. This conjecture is 
confirmed by the present calculations. 



5 Magnetoresistance effects 



Now we address the magnetoresistance property of the model Co/OSE/Co 
spin-valve system, measured as 

mr= Rap-Rp (15) 
-Rap 

Here, Rap (Rp) denotes the electric resistance with the antiparrellel (parallel) 
magnetization orientation of two FM electrodes. To study that, we turn on the 
conductance of OSE by choosing Ep = Ep + 0.8 eV (cf. Fig. 5), and calculate 
the steady-state currents [Eq. (8) with Eqs. (10)-(12)]. The resulting I-V char- 
acteristics of the model Co/OSE/Co spin- valve device in the conducting state 
for both the P- and AP-configurations of electrodes are shown in Fig. 6, and 
the corresponding MR behavior in the insert. The I-V characteristics for the 
P-configuration exhibits roughly a linear behavior (with Rp tv 17M.fl) and can 
be considered as a metallic phase. In contrast, for the AP-configuration the 
current increases much slower and is nonlinear as the bias voltage increases. 
The resulting MR shown in the insert of Fig. 6 decays nearly exponentially as 
the applied bias voltage increases from V = where MR« 100% to V = 0.5 V 
where MR = 60%. The relatively large MR here in comparing with the typical 
experimental values may partially be attributed to the fact that our model 
system does not include spin-flip processes. As shown in Fig. 6, -Rp is almost 
constant for the range of voltage reported here; thus, the bias-dependence of 
MR is governed by the .Rap in the AP-configuration, where the charge carri- 
ers in V > are the up-spin electrons driven from the full-filled (majority) 
band of the L-electrode to the half-filled (minority) band of the R-electrode 
(cf. Fig. 1). As the up-spin electrons traverse the polaron level of the OSE 
barrier, the electric resistance Rap decreases progressively as the bias voltage 
increases, resulting in the MR of an appearance of exponential decay. 

We now turn to the effect of OSE length on the transport, which has been 
considered so far only with the fixed value of N a = 10. The experimental con- 
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Fig. 6. The I-V characteristic of the OSE spintronics at the conducting state, with 
the electrochemical potential of Ep = Ep + 0.8 eV, in the parallel (solid curve) and 
antiparallel (dashed curve) configurations. The insert is the MR as the function of 
bias voltage. 

ductance often decreases exponentially as the length of the molecular chain 
increases [47]. The length-dependent conductance can be investigated theo- 
retically with the help of the equilibrium Green's function formalism [48], and 
the exponential length dependence of the conductance characterizes the co- 
herent transport [49] . Figure 7 depicts our theoretical results using the NEGF 
formalism at the bias voltage V = 0.4 V. The evaluated conductances do show 
near exponential dependences, for both the P- and the AP-configuration. The 
latter decays faster, leading to the non-uniform OSE-length dependence of 
MR, as shown in the insert of Fig. 7. 



6 Summary and concluding remarks 



In summary, we have investigated the spin-dependent charge transport in or- 
ganic spintronics devices, on the basis of the SSH model combined with the 
NEGF formalism. We have demonstrated the polaron mechanism for the con- 
ductance switching, which can be controlled either by an appropriate bias 
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Fig. 7. The length dependence of the conductance in the parallel (solid curve) and 
antiparallel (dashed curve) configurations, at the bias voltage of V = 0.4 V. The 
insert is the corresponding length dependence of the MR. 

voltage (Sec. 3) or by changing the electrode electrochemical potential (Sec. 4). 
Despite of its dependence on the electrode property, the conductance switch- 
ing is generally a property of nanosize organic semiconductors. To have a one- 
or two-step conductance switching, it is not necessary for the lead to be mag- 
netic. It is believed that the two-step conductance switching can be observed 
even in some nanostructures with nonmagnetic or paramagnetic electrodes. 
However, using magnetic electrodes does offer an additional handle on the 
conductance control (see the discussion to the end of Sec. 3); thus it has po- 
tential applications in organic electronics. We have also studied the behavior 
of magnetoresistance (Sec. 5). The calculated magnetoresistance shows an ex- 
ponential decay as the applied bias voltage increases in the V > region. The 
dependence of the magnetoresistance on the length of molecular junction is 
also addressed. 

It is worth to mention here that a decay behavior of MR has been reported 
experimentally in an LSMO/OSE(Alq 3 )/Co spin- valve device, where the MR 
has a maximal value at V ~ and decreases exponentially as the applied 
voltage increases [21]. Similar decay behavior is also observed for OSEs with 
magnetic tunnel junctions with LSMO and Co or CoFe electrodes [50,51]. 
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These experiments are for layer devices with the OSE thicknesses that may 
lead to the spin-dependent transport of diffusive. The spin diffusion can be 
incorporated phenomenologically into the present model to account for the 
relatively fast decay behavior observed experimentally. It is however beyond 
the scope of this paper and will be considered elsewhere. 

The present calculations are done for the steady or quasi-steady state electric 
transport. The dynamical polaron effects that are not considered in this work 
may be importance for the transport in organic spintronics, especially in the 
presence of a time-dependent external field. Some theoretical investigations 
have been devoted to this issue. The method used there is the wave-packet 
evolution approach [52] or non-adiabatic time-dependent Schrodinger equation 
method [53]. The resulting transport mechanism depends rather sensitively 
on the time scale of the electronic motion than that of the lattice distortion. 
Increasing the lattice vibrational frequency, the transport can change from 
electron-like to polaron-like behavior. The strong lattice fluctuation may also 
lead to the polaron transport mechanism even in the low phonon-frequency 
regime [52] . Nevertheless, the polaron mechanism studied in this work for the 
hysteretic conductance switching and giant magnetoresistance is expected to 
remain at least qualitatively valid. 
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